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Abstract 

The expulsion of magnetic flux lines from a growing superconducting core of a 
quark star has been investigated using the idea of impurity diffusion in molten 
alloys. The possibility of Mullins-Sekerka normal-superconducting interface insta- 
bility has also been studied. 
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If the matter density at the core of a neutron star exceeds a few times normal 
nuclear density (e.g. > 3n , where no = 0.17fm -3 , the normal nuclear density), a 
deconfining phase transition to quark matter may take place. As a consequence 
a normal neutron star will be converted to a hybrid star with an infinite cluster 
of quark matter core and a crust of neutron matter. If the speculation of Witten 
[1] that a flavour symmetric quark matter may be the absolute ground state at 
zero temperature and pressure is true, there is a possibility that the whole star 
will be converted to a star of strange quark matter (SQM) known as strange 
star. In normal quark matter the strange quarks are produced through flavour 
non-conserving weak processes which ultimately lead to a dynamical chemical 
equilibrium among various constituents. 

From the observed features in the spectra of pulsating accreting neutron stars 
in binary system, the strength of surface magnetic field of a neutron star is found 
to be ~ 10 12 G. At the core region of a newly born neutron star it probably reaches 
~ 10 18 G [2]. In a recent publication [3] we have shown that if the magnetic field 
intensity exceeds some critical value which is the typical strength of magnetic field 
at which the cyclotron lines begin to occur or equivalently at which the cyclotron 
quantum is of the order of or greater than the rest mass of the particle considered 
or the de Broglie wave length is of the order of or greater than the Larmor radius 
of the particle, there can not be nucleation of a single quark matter bubble in the 
metastable neutron matter. The surface as well as the curvature energies diverge 
in this case. As a consequence the new stable (quark matter) phase can not be 
thermodynamically favourable over metastable (neutron matter) phase. Therefore 
to achieve a first order deconfining transition initiated by the nucleation of quark 
droplets at the core of neutron star, we assume that the strength of magnetic field 
throughout the star is much less than the corresponding critical value. In the case 
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of electron of mass 0.5MeV, this critical field is ~ 4.4 x 10 13 G, for light quarks of 
current mass 5MeV, it is ~ 10 15 G, whereas for s-quark of current mass 150MeV, 
it is ~ 10 20 G [4-6]. 

Now for a many body fermion system, the microscopic theory of supercon- 
ductivity suggests [7] that if the interaction favours formation of pairs at low tem- 
perature, the system may undergo a phase transition to a superfluid state. This is 
expected to occur in the dense neutron matter present in neutron star [8,9]. On the 
otherhand, if the particles carry charges, the paired state will be superconducting. 
In the case of an electronic system BCS theory applied to study the superconduct- 
ing properties [7] . One electron of momentum k and spin s combines with another 
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one of momentum —k and spin — s and form a Cooper pair. The coupling is me- 
diated by the electron-phonon interaction in the lattice. In the case of SQM, the 
basic quark-quark interaction is attractive at large distance and consequently the 
BCS pairing mechanism is also applicable here. For a highly degenerate system, 
which is true in strange star, the pairing takes place near the fermi surface. The 
other condition that must be satisfied to form Cooper pair is that the temperature 
(T) of the system should be much less than the superconducting energy gap (A), 
which is a function of the interaction strength and the density of the system. This 
is the most important criterion for the occurance of superconducting transition. In 
the case of strange star, only quarks can form Cooper pairs. The electrons, whose 
density is extremely low, form highly degenerate relativistic plasma, are unlikely 
to form Cooper pairs. The kinetic energy of the electronic part dominates over 
its attractive potential energy, and as a result the corresponding superconducting 
transition temperature is extremely low and may not be achieved in a strange star. 
The relativistic theory of superfluidity and superconductivity for a fermion system 
was given by Bailin and Love [10]. Recently, Horvath et al [11] and also we have 
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studied [12] the superconductivity of quark matter using the results of ref. [10]. 
We have also studied qualitatively the magnetic properties of quark matter in ref. 
[12] for massive quarks. 

Now the quarks of same fermi energy can only combine to form Cooper pairs. 
Since the u and d current masses are equal and also their chemical potentials are 
almost identical, whereas s quark is much heavier than u and d quarks and also 
its chemical potential is different, we can have only uu, dd, ud and ss Cooper 
pairs in the system. For iso-spin 1/2 flavours, the contribution may come either 
from iso-scalar or iso-vector channels. It was shown in ref. [10] that the pairing of 
a u — d system will be favoured by iso-scalar combination rather than iso-vector 
channel. On the other hand the s — s combination is a triplet state with J p = 1 + . 
Now from the conclusions of refs. [10-12] we know that if a normal SQM system 
undergoes a superconducting phase transition, the newly produced SQM phase 
will be a type I superconductor. We have also seen that the critical magnetic 
field for such pairing is ~ 10 16 G for n ~ 2 — 3n , which is indeed much larger 
than the typical pulsar magnetic field. The corresponding critical temperature is 
10 9 — 10 10 K, this can possibily be achieved in quark star, which is expected to 
be an extremely cold object. Since the quark star magnetic field is less than the 
corresponding critical field for the destruction of superconductivity, during such 
phase transition, the magnetic flux lines from the superconducting quark sector 
of the strange star will be pushed out towards the normal crust region. Unlike 
a small type I superconducting laboratory sample placed in an external magnetic 
field (jcritical value) for which the expulsion of magnetic field takes place almost 
instantaneously, in this particular case the scenario may be completely different. 

The aim of the present note is to investigate using the idea of impurity dif- 
fusion in molten alloys, the expulsion of magnetic flux lines from growing super- 
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conducting core of a strange star. We have also studied the possibility of Mullins- 
Sekerka normal-superconducting interface instability [13] in quark matter. This 
is generally observed (i) in the case of solidification of pure molten metals at the 
solid-liquid interface, if there is a temperature gradient. The interface will always 
be stable if the temperature gradient is positive and unstable otherwise, (ii) during 
solidification of molten alloys. In alloys, the criteria for stable / unstable behaviour 
is more complicated. It is seen that, during the solidification of an alloy, there is 
a substantial change in the concentration ahead of the interface. Here solute dif- 
fusion as well as the heat flow effects must be considered simultaneously. As we 
will see, the particular problem we are going to investigate is analogous to solute 
diffusion during solidification of an alloy. 

It has been assumed that the growth of superconducting quark bubble has 
started from the centre of the star and we use the nomenclature controlled growth 
for such phenomenon. If the magnetic field strength and the temperature of the 
star are much less than their critical values, the normal SQM phase is thermody- 
namically unstable relative to the corresponding superconducting one. Then due 
to fluctuation, a droplet of superconducting quark matter bubble may be produced 
in metastable medium. If the size of this superconducting bubble is greater than 
the corresponding critical value, it will act as the nucleating centre for the growth 
of superconducting quark core. The critical radius can be obtained by minimiz- 
ing the free energy and is given by r c = lQira/Bm '' 2 [1 — (B m /Bm) 2 ], where a 
is the surface tension, which is greater than zero for a type I superconductor- 
normal interface, is the critical magnetic field. In presence of a magnetic 

(c) 

field B m < Brn , the normal to superconducting transition is first order in nature. 
As the superconducting phase grows continuously, the magnetic field lines will be 
pushed out into the normal quark matter crust. This is the usual Meissner effect, 
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observed in type I superconductor. We compare this phenomenon of magnetic flux 
expulsion from a growing superconducting SQM core with the diffusion of impu- 
rities from a molten metal. The formation of superconducting zone is compared 
with the solidification of molten metal. It is known from simple thermodynamic 
calculations that if the free energy of molten phase decreases in presence of im- 
purity atoms, then during solidification they prefer to recide in the molten phase. 
In this particular case the magnetic field lines play the role of impurity atoms, 
the normal quark matter phase plays the role of molten metal, whereas the su- 
perconducting phase can be compared with the frozen solid phase. (This idea was 
recently applied to baryon number transport during first order quark-hadron phase 
transition in the early Universe, where baryon number replaces impurity, quark 
phase replaces molten metal and hadronic matter replaces that of solid metal [14]). 
The magnetic flux lines prefer to recide in the normal phase. Then the Meissner 
effect can be restated as the solubility of magnetic flux lines in the superconducting 
phase is zero (of course, there is a finite penetration depth). 

The dynamical equation for the flux expulsion can be obtained from the sim- 
plified model of sharp normal-superconducting interface. The expulsion equation 
is given by the well known diffusion equation [15] 

^ = DV 2 B m (1) 

where B m is the magnetic field intensity and D is the diffusion coefficient, given 
by 

c 2 

D = (2) 

where a n is the electrical conductivity of the normal SQM phase, for supercon- 
ducting phase B m = 0. Following ref. [16], the electrical conductivity of SQM for 
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B m = is given by 



= 5.8 x 10 



25 




(3) 



in sec -1 , where a c is the strong coupling constant and T w = T/10 10 K, the value 
of this electrical conductivity in the case of strange quark matter is ~ 10 26 sec -1 . 
We have used this expression to get an order of magnitude estimate of electrical 
conductivity of quark matter. In actual calculation one has to evaluate a n in 
presence of B m . In that case, a n can not be a scalar quantity. In particular, for 
extremely large _B m , the components orthogonal to B m tend to zero. The quarks 
can only move along the direction of magnetic field, across the field resistivity 
becomes infinity. 

A solution of eqn.(l) with spherical symmetry (which may not be a valid 
assumption) can be obtained by Green's function technique, and is given by (for 
a very general topological structure, no analytical solution is possible) 



where u± = (r ±r')/2(Dt) 1 / 2 and Bm{r) is the magnetic field distribution within 
the star at t = 0, which is of course an entirely unknown function of radial coor- 
dinate r. To obtain an idea of magnetic field diffusion time scale (td), we assume 
Bm\r) = Bm^ = constant (in reality, this is not possible inside the star). Then 
using the other approximate result for electrical conductivity (which is valid for 
zero magnetic field case), given by eqn.(3), we have T£> = 10 — 20 yrs. With this 
constant (r), eqn.(4) 





fl ra (r,f) = Bg) 



2 r \ it J 




} (5) 



where 



erf (x) 
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is the well known error function. 

From the eqn.(5), using the approximate results for electrical conductivity, 
given by eqn.(3) (which is valid for B m = 0) one can get an estimate of time 
scale for the expulsion of magnetic lines of force and is ~ 10 — 20 yrs. Latter 
we shall see that such a time scale can also be obtained from stability analysis of 
planer normal-superconducting interface. Again the time dependence of the elec- 
trical conductivity profile for normal quark matter is not known. This is another 
uncertainty in obtaining exact solution for the diffusion equation. 

Therefore, almost nothing can be said about the growth of superconducting 
zone and the expulsion of magnetic flux lines from this region by solving eqn.(l). 
We shall now study the morphological instability of normal-superconducting in- 
terface of quark matter in the star using the idea of solute diffusion during solid- 
ification of alloys. The motion of normal-superconducting interface is extremely 
important in this case and has to be taken into consideration. Then instead of 
eqn.(l) which is valid in the rest frame, an equation expressed in a coordinate 
system which is moving with an element of the boundary layer is the correct de- 
scription of such superconducting growth, known as Directional Grwoth, and the 
equation is called Directional Growth Equation, and is given by 

dB m dB m 7-.V72 d ta\ 

-qT-v-^ = D V B m (6) 

where the motion of the interface is along the z-axis and v is the velocity of the 
front. This diffusion equation must be supplemented by the boundary conditions 
at the interface. The first boundary condition is obtained by combining Ampere's 
and Faraday's laws at the interface, and is given by 

B m v\ s =-D(VB m ).h\ s (7) 

where h is the unit vector normal to the interface directed from the normal phase 



to the superconducting phase. This is actually the continuity equation for mag- 
netic flux diffusion. The rate at which excess magnetic field lines are rejected 
from the interior of the phase is balanced by the rate at which magnetic flux 
lines diffuses ahead of the two-phase interface. This effect makes the boundary 
layer between superconducting-normal quark matter phases unstable due to excess 
magnetic field lines present on the surface of the growing superconducting bub- 
ble. Local thermodynamic equilibrium at the interface gives (Gibbs-Thompson 
condition) 

Bm |.« B£ (i _ J^-\ = B W (l - 5c) (8) 
V RB$ J 

where 5 is called capillary length with a the surface tension, c is the curvature 
= 1/R (for a spherical surface), and Bm is the thermodynamic critical field. 

To investigate the stability of superconducting-normal interface, we shall fol- 
low the original work by Mullins and Sekerka [13], and consider a steady state 
growth of superconducting core, then the time derivative in eqn.(6) will not ap- 
pear. Introducing r± = (x 2 + y 2 ) 1 ^ 2 as the transeverse coordinate at the interface, 
we have after rearranging eqn.(6) 



8 2 Id 8 2 v d 

H Ti + T^T + 



dr | r± dr± dz 2 D dz 



Bm = (9) 



Following the most common approximation which is made in the case of freezing 
of molten solid is that the solidification is occuring under steady state condi- 
tion, which in this particular case is the normal to superconducting phase tran- 
sition, and that, therefore, the concentration of magnetic flux lines and normal- 
superconducting interface morphology are independent of time. The principal 
disadvantage of this assumption is that no evolution of the interface shape can 
occur. The result of this constraint is that the solution to the basic diffusion 
problem is indeterminate and a whole range of morphologies is permissible from 
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the mathematical point of view. In order to distinguish the solution which is the 
most likely to correspond to reality, it is necessary to find some additional criteria. 
Examination of the stabilities of a slightly perturbed growth form is probably the 
most reasonable manner in which to treat this situation. In the following we shall 
investigate the morphological instability of normal-superconducting interface fol- 
lowing eqn.(9). Assuming a solution of this equation expressed as the product of 
separate functions of r± and z and setting the separation constatnt equal to zero 
and using the boundary condition given by eqn.(8), we have for an unperturbed 
boundary layer moving along z-axis 

B m = B$e-* v l D = B$e- 2 *' 1 (10) 

where / = 2D/v is the layer thickness, which is very small for small D. Mthe- 
maticaly, the thickness of this layer is infinity. For practical purpose an effective 
value / can be taken. The order of magnitude estimates or limiting values for the 
three quantities D, v and / can be obtained from the stability condition of planer 
interface, which will be discussed latter. 

Due to excess magnetic flux lines at the interface, the form of the planer 
normal-superconducting interface described by the equation z = is assumed to 
be changed by a small perturbation represented by the simple sine function 

z = esm(k.r±) (11) 

— * 

where e is very small amplitude and k is the wave vector of the perturbation. Then 
the perturbed solution of the magnetic field distribution near the interface can be 
written as 

B m = B^e~ vz / D + At sm(k.f ± )e- bz (12) 

where A and b are two unknown constants. Since the solution should satisfy the 
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diffusion equation (9), we have 

b 



2D 



+ 



1/2 



(13) 



To evaluate A, we utilise the assumption that e and esm(k.f±) are small enough 
so that we can keep only the linear terms in the expansion of exponentials present 
in eqn.(12). Then after some straight forward algebraic manipulation, we have 



A - - — S (s) 
A- D B m 



(14) 



The expression describing the magnetic field distribution ahead of the slightly 
perturbed interface then reduces to 



B = B^ 

J ->m J -> m 



z/D _ 



D 



e sin(k.f±)e~ 



(15) 



Now from the other boundary condition (eqn.(8)) we have 



= B^ - 



4ttoB^ 

r(*)2 ' 
-Dm 



(16) 



where c = z /(l + z' 2 f/ 2 is the curvature at z = esin(/c.rj_) and prime indicates 
derivative with respect to rj_. 

Neglecting z' 2 , which is small for small perturbation, we have 



Btf=B£)+rk 2 S 



(17) 



where T = AnuB^ / B^ 2 and we have replaced esin(/c.rjj by S. Since the ampli- 
tude of perturbation e is extremely small, the quantity S is also negligibly small. 
Now the eqn.(17) is also given by 



B^ = B^ + GS 



where 



G 



dB w 

dz 



\z=S z 



V 

D 



1 - V -jpj fl&> - bAS(l - bS) 



(18) 



(19) 
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Combining these two eqns., we have 

k 2 r + ^{i- v -§) ^ - b ^s(i -bS) = o (20) 

This expression determines the form (values of k) which the perturbed interface 
must assume in order to satisfy all of the conditions of the problem. To analyse the 
behaviour of the roots we replace right hand side of eqn.(20) by some parameter 
—P. (We have taken — P in order to draw a close analogy with the method given 
in ref.[13]). Rearranging eqn.(20), we have 

-^r-£(i-^) B <:> + ^d-«) = p (2i) 

(in ref. [13] the parameter P is related to the time derivative e of the amplitude 
of small perturbation). If the parameter P is positive for any value of k, the 
distortion of the interface will increase, whereas, if it is negative for all values of 
k, the perturbation will disappear and the interface will be stable. In order to 
derive a stability criterion, one only needs to know whether eqn.(21) has roots for 
positive values of k. If it has no roots, then the interface is stable because the 
P — k curve never rises above the positive /c-axis and P is therefore negative for all 
wavelengths. We have used Decarte's theorem to check how many positive roots 
are there. It is more convenient to express k in terms of b and then replacing b by 
u> + v/D, which gives 

^ (r + - u (r + *f* - £ - (i - ^ - , 

(22) 

This is a quadratic equation for u>. The first and the third terms are always 
negative. The second term will also be negative if 

r+*^^-5«S>0 (23) 
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Then it follows from Decart's rule that if the condition (23) is satisfied, there 
can not be any positive root. Which implies that the small perturbation of the 
interface will disappear. Since the amplitude of perturbation is assumed to be 
extremely small, the quantity S = esin(/c.fj_) is also negligibly small. Under such 
circumstances the middle term of eqn.(23) is much smaller than rest of the terms. 
The Decart's rule given by the condition (21) can be re-written as 

r > flWfi (24) 

Which after some simplification gives the stability criterion for the plane unper- 
turbed interface, given by 

r>(s)3 q 

a > ^-T (25) 

From the stability criterion, it follows that the normal-superconducting in- 
terface energy /area of quark matter has a lower bound, which depends on the 
interface magnetic field strength, critical field strength and also on the pertur- 
bation term S. An order of magnitude of this lower limit can be obtained by 
asuming = 10~ 3 B^ . (Since the critical field Bm ~ 10 16 G, and the neutron 
star magnetic field strength B m ~ 10 13 G, we may use this equality). Then the 
lower limit is given by 



a L « 10" 9 MeV/fm 2 (A) (26a) 



(s) (c) 

On the other hand for Bm = O.lSrn , we have 



a L « 10" 3 MeV/fm 2 ( ^- ) (266) 
The approximate general expression for the lower limit is given by 

a L w h 3 MeV/fm 2 (J^ (26c) 
13 



where h = B$/B$. There for the maximum value of this lower limit is 



a£ ax - w 1 MeV/fm 2 (26d) 

when the two phase are in thermodynamic equilibrium. Of coures for such a 
strong magnetic field, as we have seen (see ref.[3]) that there can not be first order 
quark-hadron phase transition, and it should be some higher order one. 

On the other hand if we do not have control on the interface energy, which 
can in principle be obtained from Landau-Ginzberg model, we can re-write the 
stability criteria interms of interface concentration of magnetic field strength Bm , 



S 



and is given by 



m 



(27) 



S{1-%S) 

This is more realistic than the condition imposed on the surface tension a. Now for 
a type I superconductor, the surface tension a > 0, which implies 1 — 2vS/D > 0. 
Therefore, we have 2vS/D < 1, and for the typical value a n ~ 10 26 sec -1 for the 
electrical conductivity of normal quark matter with zero magnetic field, the profile 
velocity v < D/2S ~ 10~ 6 /S cm/sec ~ 1 cm/sec for S ~ 10 -6 cm. Therefore 
the interface velocity < 1 cm/sec for such typical values of a n and S to make the 
planer interface stable under small perturbation. Now the thickness of the layer 
at the interface is / = 2D/v > 10 -6 for such values of D (or a n ) and v. Here 
S is always greater than 0, otherwise, the magnetic field strength at the normal- 
superconductor interface becomes unphysical. As before, if the second term of 
eqn.(23) is negligibly small compared to other two terms, we have 



B (s) 



III. ^ 



1/3 



(28) 



Therefore we may conclude by saying that if a superconducting tnasition can 
take place in a quark star, the magnetic properties of such bulk object are entirely 
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different from that of a small laboratory superconducting sample. Expulsion of 
magnetic flux lines from the superconducting zone is not instantaneous. The 
typical time scale is 10 — 15 yrs. Due to the presence of excess magnetic flux lines 
at the interface, which is possibly true if the diffusion rate of magnetic lines of 
forces in the normal phase is less than the rate of growth of the superconducting 
zone, the characteristic of normal-superconducting boundary layer may change 
significantly. Of course, it depends on the magnitude of surface tension a. It may 
take dendritic shape instead of planer structure. The stability of planer interface 
also depends on the strength of interface magnetic field, if we do not have control 
on the interface energy and are given by eqns.(27) and (28). How to get an 
experimental evidence for such an unusual shape is a matter of further study. 



15 



References 

1. E. Witten, Phys. Rev. D30, 272 (1984). 

2. J. Trumper et al, Ap. J. 219, L105 (1978); W. A. Wheaton et al, Nature 272, 
240 (1979); D. E. Gruber et al, Ap. J. 240, L127 (1980); T. Mihara et al, 
Nature 346, 250 (1990). 

3. S. Chakrabarty, Phys. Rev. D51, 4591 (1995). 

4. S. Chakrabarty, Astrphys. and Space Sci., 213, 121 (1994). 

5. S. Chakrabarty and A. Goyal, Mod. Phys. Lett. A9, 3611 (1994). 

6. S. Chakrabarty and P. K. Sahu, Phys. Rev. D (in press 1996); S. Chakrabarty, 
Phys. Rev. D (submitted). 

7. A.L. Fetter and J. D. Walecha, Quantum Theory of Many Particle System, 
McGraw Hill Book Company, New York, 1971. 

8. G. Baym, C. Pethick and D. Pines, Nature 224, 224 (1969). 

9. M. Baldo et al, Nucl. Phys. A536, 349 (1992). 

10. D. Bailin and A. Love, Phys. Rep. 107, 325 (1984). 

11. J. E. Horvath et al, Mod. Phys. Lett. A7, 995 (1992). 

12. S. Chakrabarty, Can. J. Phys. 71, 488 (1993). 

13. W.W. Mullins and R. F. Sekerka, Jour. Appl. Phys. 34, 323 (1963); 35, 444 
(1964). 

14. F. C. Adams, K. Freese and J. S. Langer, Phys. Rev. D47, 4303 (1993); Marc 
Kamionkowski and K. Freese, Phys. Rev. Lett. 69, 2743 (1992). 

15. J. S. Langer, Rev. Mod. Phys. 52, 1 (1980). 

16. P. Haensel and A. J. Jerzak, Acta Phys. Pol. B20, 141 (1989); P. Haensel, 
Nucl. Phys. B24 (proc. of the Int. Workshop on Strange quark Matter in 
Physics and Astrophysics, Univ. of Aarhus, Denmark, May 20-24, 1991), 23 
(1991). 



16 



